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Abstract 

We study diffractive pliotoproduction of J /ip by taking the charm quark as a heavy quark. 
A description of nonperturbative effect related to J/ip can be made by using NRQCD. In 
the forward region of the kinematics, the interaction between the cc-pair and the initial 
hadron is due to exchange of soft gluons. The effect of the exchange can be studied by using 
the expansion in the inverse of the quark mass rric- At the leading order we find that the 
nonperturbative effect related to the initial hadron is represented by a matrix element of field 
strength operators, which are separated in the moving direction of J/ip in the space-time. 
The S-matrix element is then obtained without using perturbative QCD and the results are 
not based on any model. Corrections to the results can be systematically added. Keeping 
the dominant contribution of the S-matrix element in the large energy limit we find that 
the imaginary part of the S-matrix element is related to the gluon distribution for x ^ {) 
with a reasonable assumption, the real part can be obtained with another approximation or 
with dispersion relation. Our approach is different than previous approaches and also our 
results are different than those in these approaches. The differences are discussed in detail. 
A comparison with experiment is also made and a qualitative agreement is found. 

PACS numbers: 12.38.-t, 12.39.Hg, 13.60.-r, 13.60.Le 
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1. Introduction 



It is usually believed that the nonperturbative QCD plays an important role in diffractive 
processes and one cannot use perturbative QCD to describe them. Recently it was pointed out[|I| 
that for diffractive production of a vector meson V like 

Y + h^h + V (1) 

can be handled with perturbative QCD provided that the virtuality of the initial photon is 
large. This enables us to make testable predictions for the process and it provides an interesting 
way to study the nonperturbative nature of the initial hadron, e.g., the structure function of 
hadrons and of nuclei. For V = J/ip it is also studied in p with perturbative QCD. Theoretically 
it is proved that the S-matrix element can be factorized[^]. Neglecting higher orders of the 
S-matrix element consists of the light-cone wave function of V, skewed parton distributions of h 
and a hard scattering kernel, the hard scattering kernel can be safely calculated with perturbative 
QCD and it is free from infrared singularities. It should be noted that the factorization also holds 
if one replaces the vector meson V with a spin-0 meson, or with a photon, the so called deeply 
virtual Compton scatter ing[|]. 

In this work we study the diffractive photoproduction of J/ip 

-f + h-^h + J/ij. (2) 

Neglecting a possible cc-content of h, the process can be imagined as the following: The photon 
splits into a cc-pair, after interactions with the hadron h through gluon exchanges the cc-pair is 
formed into J/ip. Because the initial photon is real, i.e., = 0, the factorization proved in 
does not apply here. If the total energy is sufficiently large, the exchanged gluons are soft and 
they can not be handled with perturbative QCD. But the charm quark can be taken as a heavy 
quark, for emissions of soft gluons by heavy quarks the heavy quark effective theory(HQET) can 
be used0, a systematic expansion in the inverse of the charm quark mass rric can be employed 
to study emissions of soft gluons. Taking the charm quark as a heavy quark it also allows us to 
use nonrelativistic QCD(NRQCD) to describe nonperturbative properties of J/ip- As an approxi- 
mation one can take J/ijj as a bound system of a c- and c-quark, in which the c- and c-quark has 
a momentum which is half of the momentum of J/ip- Corrections to this approximation can be 
systematically added in the framework of NRQCD0. 

The process studied here has a close similarity to the decay of J/ip 

J/?/' — i> 7* + TT + TT — > + + TT -|- TT (3) 

in the kinematic region where the pions are soft. This decay is studied in J^, §. The exchange 
of soft gluons between the cc pair in J/tp and the pion pair is responsible for the decay. In ^ 
it is shown that one can use a technique of path integral for the exchange of soft gluons without 
invoking perturbative QCD. It is also shown that at the leading order of the decay amplitude 
derived with the technique of path integral is the same as that derived by assuming that the 
exchange is of two soft gluons in a special gauge. In this work we will take a suitable gauge and 
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assume the two-gluon exchange in the gauge to derive the S-matrix element for the diffractive 
process. The assumption can be justified by using HQET: In a suitable gauge the probability for 
a c-quark emitting 1 or 2 gluons is proportional to m~^, while the probability for emission of more 
than 2 gluons is at order of m~" with n > 1. It is interesting to note our result can be derived 
without taking a suitable gauge and the assumption of two gluon exchange. This can be done 
with the technique of path integral as that derived for the decay in we briefly sketch how to 
derive our result in this way and details may be found in [0]. Because the result is derived without 
the assumption of two gluon exchange in a suitable gauge, our result actually includes effects of 
exchange of more than 2 gluons in an arbitrary gauge. This will be discussed in detail after our 
result is represented. The technique of path integral has been used to study exchanges of soft 
gluons between two light quarks with large momenta[§], where the exchanges were responsible for 
diffractive scattering of light hadrons. 

Our results for the S-matrix element consists of a NRQCD matrix element and a distribution 
amplitude of gluons in the light hadron h. The NRQCD matrix element represents the nonper- 
turbative effect related to J/ip, the distribution amplitude is deflned by a matrix element of two 
field strength operators separated in the moving direction of J/ip in the space-time. It should be 
emphasized that the obtained results are not based on any model, corrections to the results can be 
systematically added in the framework of QCD. At the leading order we consider, the produced 
J/ip has the same polarization of the photon, i.e., the produced J/ip is transversely polarized. 
In the limit of large beam energies, i.e., s oo, the dominant contribution of the amplitude is 
related to the skewed gluon distribution of h. With a reasonable assumption the forward S-matrix 
element can be related to the usual gluon distribution gh{x) for a; 0. This enables us to pre- 
dict the forward differential cross section with available information of gh{x)- With this result it 
provides an interesting way in experiment to access the small x-region of gh{x)- 

There exist two approaches for the diffractive process in Eq.(l). One is to use perturbative 
QCD[||, 0, with several approximations one obtains the S-matrix element related to the gluon 
distribution. Another one is based on the perturbative QCD result for interaction of a small 
transverse-size dipole of a quark pair, which is formed into V. The interaction of the dipole with 



the initial hadron is through two-gluon exchange [pl]|] . Both approaches have close similarities. In 



these approaches all hadrons including J/ip should be taken as massless at the leading order of 
Q^^. Because an expansion in is used, one can not take the limit ^ to obtain predictions 
for the photoproduction. The expansion in also implies that the S-matrix element is obtained 
in the limit of s — > oo since Q'^/s remains finite. In |[Tl[] the approach of the dipole interaction |10 



is used, and the effect of the nonzero mass of J/ijj is taken into account. It is found that one can 
take — to have the S-matrix element for the photoproduction, where the S-matrix element 
is also related to the gluon distribution. Similar results are also obtained in p| . It is questionable 
if the limit ^ can be taken because the higher orders oi Q are neglected. Our approach is 
distinctly different. We start directly from the process in Eq.(2) and obtain the S-matrix element 
for a moderate s. Then we take the limit s oo, and the forward S-matrix element in the limit 



is related to the usual gluon distribution. Our results are also different than those given in g, |TI 
We will discuss the differences in detail. 

At first look, one may generalize our results to the case where the initial photon is virtual 
with a small Q^. For transversely polarized J/ip the generalization is straightforward. But, for 
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longitudinally polarized J/V' the generalization seems not possible, because the quark mass rric is 
involved in the polarization vector of J/ip, this can spoil the expansion in m~^. The production 
of longitudinally polarized J/ ip deserves therefore a further study and we will briefly discuss the 
problem. 

Our work is organized as the following: In Sec. 2. we introduce our notations and derive the 
S-matrix element in the diffr active region at the leading order of mj^. The result is derived by 
taking the exchange of two soft gluons in a special gauge into account. Wc will briefly discuss how 
to derive it with the technique of path integral. We will also briefly discuss the problem of the 
U{1) gauge invariance and of the production of longitudinally polarized J/ip. In Sec. 3. we derive 
the forward S-matrix element in the large energy limit. The S-matrix element is related to the 
usual gluon distribution Qhix) with x — > 0. We discuss the difference between our approach and 
another approach. In Sec. 4. we compare our results with experiment. Sec. 5 is our summary. 



2. The soft gluon approach 



We consider the process 

7(A;) + h{p) ^ h{p + A) + J/tlj{k - A), (4) 

where the momenta are given in the brackets. The Mandelstam variables are defined as 

s = (A;+p)2, i^A^. (5) 

We study the process in the kinematic region where \t\ is at order of A.qqj-) and each component 
of A is at order of Aqcd- h is any light hadron whose mass m is at order of Aqcd- By taking the 
charm quark as a heavy quark we have 

Mj/^>|t|, s>|i|. (6) 

It should be noted that we do not require that s ^ ^j/^j our result presented in this section can 
be applied to a wide range of s. But the value of s should be not too small so that \t\ can be small 
enough. The smallest value of |i| can be approximated by: 

hence the value of s should satisfy: 

s{s - Mj/^) » m^Mj/^. (8) 

It should be noted that at the threshold we can also have the conditions in Eq.(6), but the 
exchanged gluons are not soft, hence, our approach presented in this work cannot be used for 
production at the threshold. 

Although our result can be derived exactly by using the technique of path integral as explained 
in the introduction, we derive here our result by assuming two-gluon exchange in a special gauge. 
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Figure 1: One of six Feynman diagrams for elastic photoproduction of J/ip. 



because the derivation in this way is straightforward. We will briefly discuss how to derive our 
result without the assumption. The contributions of two-gluon exchange to the S-matrix element 
can be represented by diagrams, one of them is given in Fig.l. The S-matrix element with two- 
gluon exchange can be obtained directly: 



{f\S\i) = -J d^xd^yd^xid'^X2 J 



d'^qi d'^q2 d^ki d*k2 
(2^(2^(2^(2^ 



■£,(fc)Ag''''^^(fci, k2. gi, g2)(27r) V(A; + fc^ - ^2 - gi - ga) 

■e-^''-^-*^-^(W|c.(x)c,(y)|0)-e^'=-^-*^-^^(/i(p + A)|G»(xi)G^(a:2)|/^(p)), (9) 



where e{k) is the polarization vector of the photon, c{x) is the Dirac field of the c-quark, the index 
i and j stands for color- and Dirac indices. A'^j''''^'^{ki,k2,qi,q2) is the scattering amplitude for 
the process: 

7(A;) + G*iku a) ^ G*(A;2, b) + c*(gi) + c*(g2), (10) 

where quarks and gluons are not necessarily on-shell. The matrix {J/il)\ci{x)cj{y)\0) represents the 
nonperturbative effect related to J/ip. It should be noted that for exchange of arbitrary numbers 
of gluons the same matrix appears in the S-matrix element. Because we take the charm quark as a 
heavy quark, the c- or c-quark in J/tp carries roughly the half of the momentum of J/ip, the effect 
induced by the deviation from the half momentum is suppressed, the suppression parameter is the 
velocity Vc of the c- or c-quark in J/ ip in its rest frame. This fact can be realized by boosting the 
moving frame of J/ip to its rest frame, in the rest frame one can then uses NRQCD to perform an 
expansion in fc[§]. We will treat the matrix in the moving frame by using HQET, and then the 
nonperturbative effect is represented by matrix elements defined in HQET. These matrix elements 
can be related to those defined in NRQCD. Although the expansion parameter in HQET and in 
NRQCD is different, but at the orders we consider here this will not cause problems. In this work 
we take nonrelativistic normalization for heavy quark states and for J/ if) state. 
We define the velocity v oi J/ip as 
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the Dirac field c{x) can be expanded in ^ with fields of HQET: 

c{x) = e-''"^^-^ + ^7 • Drhix)^ + e+''"^''-^ + ^7 • ^t^(^)} + Oimf), (12) 

where Dj. = — f ■ ^, is the covariant derivative. h{x) and are fields of HQET, h{x) 
can only annihilate a heavy quark and g{x) can only create an heavy antiquark. These fields have 
the property 

■vh{x) ^h{x), ■vg{x)^ -g{x), (13) 
they also depend on the velocity v. With these fields the matrix can be written 

(J/^|Q(x)c,(y)|0) = e~(-+^)(J/^|/z,(x)^,(y)|0) + • • • (14) 

where the • • • stand for higher orders in mj^. With the expansion the c- and c quark in J/%1) 
carries the momentum rricV plus some residual momentum, the effect of the residual momentum 
is represented by the space-time dependence of the matrix element of HQET fields. Because the 
effect is small, we can neglect the dependence. We obtain: 
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imcV-{x+y) 

6' 



(1-^^.^.^^^(1+^ 



{J/,p\h^-e{v)g\Q) 



+ ..-, (15) 

where the matrix is diagonal in the color-space, the matrix element is of local fields in HQET, 
e*{v) is the polarization vector of J/ip. We will neglect the higher orders represented by • • • in 
the above equations, then the momentum of J/'^ is approximated as 2mcV. The matrix element 
{J / ip\h{x)^ • e{y)g\Q) is related to a NRQCD matrix element defined in the rest frame of J/V'- The 
relation reads: 

\/^{jmh-i ■ e{v)g\Q) = -(J/^I^V ■ £(v = 0)x|0) (16) 

where ip and x ^I'c NRQCD fields for the c- and c quark respectively, (jj (i = 1, 2, 3) is the Pauli 
matrix. This NRQCD matrix element can be determined from the leptonic decay of J/ip 

r(W - e+e-) = aLQ^^KWlV'V • s{v = 0)xm'. (17) 

Through examination of contributions of higher orders in Eq.(15) and Eq.(14) one may find that 
they are suppressed by v"^ relatively to the leading order contribution. 
With the result in Eq.(15) the S-matrix element reads: 

{f\S\t) = \ j d'x.d'x, j ^^(27r)^5^(fc + 2m,v) 

' ^ :(J/V^|V^to- . £(v = 0)x|0)) • {~eQ,g%,R''''{m,v, h, k^)) 



" 6V^^ 2 
.e^k^-^^-^k2-^^{h{p + ^)\Gl{x)Gl{y)\h{p)), (18) 
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where 

- \eQ,gl5abR^''{m,v, h, h) = Tr i^e,{k)A^p"''' {h, h, m,v, nicv)^ ■ ^*(^)^^^} • (19) 
Now we take the special gauge 

V ■ G{x) = 0. (20) 

In the gauge we have: 

v,G>^'^{x)^v,^G%x), (21) 

where G^^'^{x) is the field strength tensor of gluon. In R^^{mcV, ki, ki and k2 are the momenta 
carried by the two exchanged gluons, their components arc small in comparison with rric, because 
the two gluons are soft gluons. Hence R'^'^{mcV, ki, k2) can be expanded in m~^, a formal expansion 
in leads to: 

rric V ■ [ki — k2) [v ■ ki — iO+){v ■ k2 + tO+) 

+ (terms proportional to or v^), (22) 

where 0^ denotes an infinitesimal positive number, which comes from quark propagators. The 
leading order is m~^. Substituting the result of R'^^{mcV,ki,k2) at the leading order in the S- 
matrix element in Eq.(18), the integrations of the components of k\, k2, xi and X2, which are 
transverse to v, can be easily performed, and the result for the S-matrix element at the leading 
order can be obtained straightforwardly. To present the result we define 

e-'^^^-^'v.v^ihip + A)\G^'^''iTv)G-'''^i-Tv)\hip)), 

-oo 2,71 

/oo rjj- 
—e'^'"'^v^v^{h{p + A)\G''^i'P{Tv)G'''''J-Tv)\h{p)) (23) 
-oo 2tt 

where we have used 

vA^-m,(^l + ^ j ^ -m,. (24) 

The variable z is related to the momentum ki and k2 by 

V ■ ki^ ~v ■ A{1 + z), v-k2 = -v-A{l-z). (25) 

The function Fr^z) are zero for \z\ > Zq with zq = {2v ■ p + v ■ A) / rric because of the conservation 
of momentum. With the function the result for the ^S-matrix element reads: 

{f\S\t) = ^ea(27r)V(fc-A-2me^)"^^^''*^'^^ 



{Jli^Wa- ■ e(v = 0)x|0) [ I 

{vAf J (l + ^-iO+)(l-^-iO+) ^ 

2^ .^.O.N4.4/,. A ..Ak)-^*i^) 



,2 



eQc{2T:Y5\k - A - 2mcv) - 



•U/^I^V . .(V = 0)x|0> (26) 



The above results are derived with the assumption of two-gluon exchange in the gauge v ■ G = 0. 
In this gauge the polarization vectors of the two exchanged gluons are perpendicular to v. To 
maintain the color-gauge invariance in other gauges a gauge link must be supplied between the 
field strength operators in Eq.(23), with the gauge link the effect of exchanges of gluons, whose 
polarization vectors are proportional to v and whose number is unlimited, is also included. Our 
result derived in this way may be unsatisfied, because the assumption of two-gluon emission sounds 
that we performed an expansion in Qs for soft-gluons and we add the gauge link by hand. It is 
possible that the cc pair emits soft gluons whose polarizations are all proportional to v and this 
emission is not suppressed by m~^. This type of contributions was excluded with the gauge. 
However our result can be derived in an arbitrary gauge without the assumption of two-gluon 
exchange. The derivation is similar to this for the decay in Eq.(3), we will briefly describe the 
derivation here, details can be found in [|^. 

Considering the process with exchange of arbitrary number of gluons in an arbitrary gauge, 
the matrix element in l.h.s. of Eq.(15) always appears in the S-matrix element. With the approx- 
imation in Eq.(15), it is equivalent to consider photoproduction of a cc-pair, where the c- and c 
quark is on-shell and has the same momentum rricV. Using the standard SLZ reduction formula 
we related the S'-matrix element to Green's functions, which can be calculated with QCD path 
integral. Imaging that we perform first the integration over c-quark fields, then the problem is 
formulated as to solve the wave functions of c- and c quark in the process under a background of 
gluon fields: 

7 ^ c + c. (27) 

The background fields vary slowly with the space-time, reflecting the fact that the exchanged 
gluons are soft. The wave functions can be solved with an expansion in m~^. At leading order, i.e., 
at order of m°, the wave- functions are obtained by multiplying the wave-functions in the free case 
with gauge links determined hj v -G. This means that at the order of m° only those gluons whose 
polarization is proportional to v, are exchanged. Because of the symmetry of charge conjugation 
these gauge links do not lead to any physical effect, i.e., the S-matrix element is zero at m°. Solving 
the wave-functions at order of one obtains exactly the same results given in Eq.(26), and the 
gauge link, which needs to be added by hand in Eq.(23), is automatically generated. Because the 
results are derived without the assumption of two-gluon exchange, they are nonperturbative. The 
results indicate that the exchange consists of two gluons, whose polarizations are transverse to the 
moving direction of J/ip, and of any number of gluons, whose polarizations are proportional to v. 

Our results show that the produced J/ip is transversally polarized at the order we consider, 
the production of longitudinally polarized J/ip is suppressed. As they stand, the results do not 
respect to the gauge invariance of electromagnetism. This can be seen that the S-matrix element 
in Eq.(26) is not zero if e{k) is replaced by k. To study the problem, we note that the performed 
expansion in is a formal expansion, the true expansion parameter is {rricV-ki)^^ or {mcV-k2)~^, 
this can also be realized by inspecting the HQET lagrangian. To identify the expansion parameter 
more clearly, we scale the momenta: 

k, = Xh, k2 = A^2, (28) 

where the components of ki or k2 are C(l), A is proportional to Aqcd and is small. We note that 
the f-dependence in R^'^ appears though v-ki and v-k2, if we use k-e{k) = {2mcV+k2 — ki)-e{k) = 0. 
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We scale these factors as: 

V ■ ki — rricV ■ ki, v ■ k2 — rricV ■ k2, 
where v ■ ki and v ■ ^2 are 0{1). Now we can expand R'^" in A, the result reads: 



(29) 



R'^''{mcV, ki, ^2) = Sp^^r" + C'( — ) + (terms proportional to or v") 



(30) 



'C 



The first term is identical to the term in Eq.(22), which is at order of A*^, the second term has 
a length form and is at order of A, this implies that the corrections to our results in Eq.(26) is 
at order of C'(;^), where corrections come not only from the term at order of C(^), but also 
from emission of more than 2 gluons in the gauge. With this examination the gauge invariance 
is violated at order of (9(^), i.e., at the next-to-leading order which we neglect, by noting that 
k ■ s*{v) = 0{\). To restore the gauge invariance, one needs to analyze the contribution at the 
next-to-leading order. Retaining only the leading order, the gauge invariance holds. 

Our results show that the produced J/i/j has the same helicity of the initial photon, i.e., 
the produced J/t/j is transversally polarized. This is easy to be understood by noting that the 
exchanged gluons at the considered order do not change the helicity of the c- or c quark. We can 
formulate our results as: 



where A^ and Xj is the helicity of the photon and of J/ip, respectively. This is our main result 
of this section, where the order of errors is also given. The first error is from neglecting higher 
orders in the mj^ expansion for emission of soft gluons, while the second is from neglecting the 
relativistic correction related to J/^. 

If the initial photon is virtual and the virtuality is small, the exchanged gluons arc also 
soft. In this case our approach can be used. It is straightforward to generalize our results to 
the production of transversally polarized J/tp with the transversally polarized photon. But, the 
generalization may not be done for longitudinally polarized J/ip with longitudinally polarized 
photon. The reason may be seen from the expansion in for i?^*" in Eq.(22). If the photon and 
J/i/j are transversally polarized, their polarization vectors have components which are all at order 
of 0{1). Then tjIc is the only large parameter in and an expansion in can be performed. 
This fact is used in the expansion in Eq.(30). If they are longitudinal polarized, their polarization 
vectors can have components which are very large in comparison with 0{1), this may prevents us 
from an expansion in mj^ for R^". It deserves a further study of the production of longitudinally 
polarized J/ip. 

3. The forwcird S-matrix element in the hmit of s ^ 00 



{f\S\z) = {2nyS\k-A-2m,v)^,^,^■{-^^eQ,){J/iJ\iJ^(T■e{^r = 0)x\0) 




(31) 
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In this section we discuss the S-matrix element in the hmit of s oo. In this limit, the 
function Fr{z) can be directly related to the skewed gluon distribution, whose definition can be 
found in 0, This can be realized by that in the limit the dominant part of v is proportional 
to a light cone vector defined below and the correlator defined in Eq.(23) can be expanded with 
operators classified with twist. The leading order is determined by twist-2 operators. However, 
the skewed gluon distribution function is not well known and this prevents us from numerical 
predictions. But, as we will see, there is a possibility to relate the imaginary part of the forward 
S-matrix element to the usual gluon distribution under certain approximation as in the case of the 
process in Eq.(l)|l|]. We will show that the imaginary part can be related to the gluon distribution 
with a reasonable assumption and the real part can be estimated by an approximation, then we 
obtain the forward S-matrix element determined by usual gluon distribution function gh{x) of h 
for a; ^ 0. Before taking the limit s oo, we note that the integral over z in Tr in Eq.(31) can 
be performed analytically. Because Fji{z) = for \z\ > Zq, we can extend the integration of z and 
exchange the integration of z and that of r in Fr{z). Using an contour integration we obtain: 

/ + / ^^^'^^^(P + A)|G"^'^''(rt;)G'^';(-rt;)|Mp)) 

JO ZTT J-oo ZTT I ^ 

CO 



dr 

/o 2^^ 

For convenience we take a coordinate system in which the photon moves in the — 2;-direction 
and the hadron h in the 2;-direction. We introduce a light-cone coordinate system, components of 
a vector A in this coordinate system are related to those in the usual coordinate system as 



(33) 



The momenta in the process can be approximated in the limit s ~ 2/c p+ oo as: 
F = (0,fc-,OT), 

2 

= ,Ot)~(p+,0,Ot), 
/ 1 - 

Af" = 

where v is approximated as a light cone vector. Using these approximated momenta and ^j/^ 
Ami > 1^1 Tr can be written: 




Tn^{2^^)gl^ T ^e-^'-^"(% + A)|G'^'+''(^n)G'^'%(-\)|Mp)), (35) 
4mc Jo Iir 2^2 

with 



2m^ 

72^^ = (0,1,Ot), x, = ^, (36) 
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where Tr is approximated by the matrix element of the twist-2 operator and corrections can be 
parameterized with higher-twist operators and they are suppressed by large scales like s. We 
consider the forward case, i.e., t —>■ 0. For t ^ and A'^ the integral in Eq.(35) can be 
approximated by the replacement: 

+ A)|G'^'+''(^n)G%(-^n)|Mp)) ^ = /(A), 

k~ r°° dX 



' Arric Jo 2t[ 

with an assumption which we will discuss later in detail. In Eq.(37) we introduce the notation 
/(A) for the forward matrix element, it has the property /(A) = /(—A). Again, /(A) also appears 
in the definition of the gluon distribution in h, in the light cone gauge we use the definition 



is 13 



^ I -i\p+x . 



2 fix 

= --/ i^cos(Ap+a;)/(A). (38) 
p+ JO ivr 

Therefore the forward S-matrix element can be related to the usual gluon distribution. At first 
look, it seems that there is an ambiguity in relating Tr with the usual gluon distribution. One 
can use the translational covariance to shift the variable A in Eq. (35), 

(Mp + A)|G»'+''(\)G%(-^n)|Mp)) 
= e-*^"^^'^(/i(p + A)|G"'+^(0)G"'+^(-An)|/i(p)). 
If one makes the approximation, 

{hip + A)\G^'+^iO)G'^'\i-Xn)\hip)) ^ {h{p)\G'^'+^iO)G'^'\i~Xn)\h{p)) = /(A), 
then one will obtain: 



{27,r)gl-— -le-^^'^^P^f{X), (39) 



dX 

ii^[27n)g;- / — 

Arric JO 2vr 

i.e., one will get a different result, because Xc in Eq.(37) is replaced by 2xc in Eq.(39) now. This 
will result in that the forward S-matrix element will be related to the gluon distribution gh{x) at 
X = 2xc instead of gh{x) at x = Xc, as we will see below in Eq.(41) and Eq.(51). Then we will 
have two different results of one theory. In the following, we will use the approximation Eq. (37) 
to derive our results and will demonstrate later that the approximation leading to Eq.(39) is not 
consistent because some neglected contributions in this approximation are at the same order as 
those kept in the approximation. 

To study the relation of Tr to the gluon distribution in detail, we write Tr as a sum of the 
integrals: 

TnU^o = W.^ if-dX^r-d_X_rOdX) (40) 
Arrir J-oo 2% Jo 2n J-oo 2n 
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where the first term is simply proportional to the gluon distribution with x — Xc-, the second and 
third terms can be combined into one integral by using the property /(A) = /(—A): 

^Rlt^o = ~^^9l^ [x9h{x)] \x=xc + '^ali^ ^ sin(Ap+Xc)/(A). (41) 

In the above equation the first term is the imaginary part of the S- matrix element, the second 
term is the real part and it can be estimated with dispersion relations. Here we estimate it by 
another method. Wc note that the limit s — > oo implies Xc 0, the asymptotic behavior of gh{x) 
with X — > is expected to be 

xgtix) ~ a;-", (42) 

where a > 0. With this behavior one can expect that the dominant contribution to the second 
term is determined by the asymptotic behavior. It is clearly that this behavior is determined by 
the behavior of /(A) for A — > oo. For A ^ oo the function /(A) goes to zero, if /(A) converges to 
zero fast enough with A — > oo, the singular behavior in Eq.(42) will not appear. If we calculate 
/(A) with perturbative theory, the result looks hke 



/(A)| 



pert. 



A"'(co+Ec„(ln(A/.))"j+---, (43) 



where jj, is the renormalization scale. The terms represented by ■ ■ • will result in singular terms, 
like (5-functions in the perturbatively calculated gluon distribution. These terms are irrelevant in 
our case. Neglecting the terms with ln(A/x), /(A) behaves like the power behavior A~^. However, 
if one sums all terms with the logarithm ln(A/i), the power behavior will be changed. Also 
nonperturbative effects will definitely change the behavior. We assume that the function /(A) 
takes the form for A — > cxd: 

/(A) ^ fasiX) = 5 (44) 
and rewrite the integral in the definition of gh{x) as: 

xgH{x) = -- cos(Ap+x)/..(A) - - / cos(Ap+x)(/(A) - /..(A)). (45) 

In the second term the function (/(A) — /as(A)) will converge to zero faster than /as(A), when 
A — oo, hence the most singular term of xgh{x) for x ^ comes from the first term. We write 
the distribution as: 

xgh{x) = Gas{x) ■ (1 + Gr{x)) , Gas{x) = ^, (46) 

where Gr{x) with a; — > 0. It should be noted that in general Gr{x) ■ Gas{x) can be singular 
for X — > 0, Gas{x) is the most singular term in gh{x) for x ^ 0. With these notations the most 
singular term in gh{x) is related to /as (A): 

r'^ dX 



Gasix) = — - / — cos(Ap+x)/„s(A) 
p""" Jo Ztt 

- -^ixp^r^\f^X-^cos{X). (47) 
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It should be noted that the integral is finite as discussed in the appendix. Comparing with Gas (x) 
in Eq.(47) we obtain the relation between a and /3 and that between Ag and a^: 



a 



1-/3, 



(48) 



With the above discussion we can realize that the dominant contribution to Tr is obtained by 
replacing /(A) with /as (A), and this dominant contribution is determined by the most singular 
term in the gluon distribution. We will only take this dominant contribution. With the form of 
/as (A) it is straightforward to calculate Tr for t —> 0: 



k- 



dX 



Arric Jo 27r 



/as (A) 



Arric 



[XcP 



poo 

/ dAA-^(cos(A) -isin(A)). 
Jo 



With the relations in Eq.(48) and results in the appendix for the integrals we obtain: 



1 — i tan(-Q;7r) 



(49) 



(50) 



For a enough small Xc one may replace Gas{xc) with [xcgh{xc)] as a good approximation: 



1 — i tan(-Q;7r) 



(51) 



In the above result the spin of the light hadrons are the same. With the S-matrix element the 
forward differential cross section reads: 



da , 



dt Stts^ 



2 riJ/i/j ^ e+e-) 1 



ttem rri] 
271^0^ T{J/ijj e+e 
3 aemrricS^ 



1 2 

gh{xc){l - Uan(-a7r)) 



(52) 



where is the summation over spin of the final hadron h and the spin average of the initial 
hadron. This result is our main result in this section. It should be emphasized that we only used 
the asymptotic behavior of gh{x) to estimate the real part of the amplitude, i.e., the term related 
to tan(|Q;7r), and the imaginary part is determined without using the asymptotic behavior, as it 
already stands in Eq.(41). With the result in Eq.(51) and Eq.(52) a problem may arise if a — > 1. If 
a is really close to 1 or equals to 1, then the cross section will become infinitely large. However, this 
can not be the case because a = 1 implies that the second moment of gh{x), which is the average 
of the momentum fraction carried by a gluon in the hadron h, is infinitely large. Experimentally 
the extracted second moment is smaller than 1. The value of a from different parameterization 
of gh , relevant to our case, is found to be 0.15 ~ 0.35 for J/^, the corresponding value for T 
production is 0.30 ~ 0.48. 
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Now we are in the position to discuss the assumption leading to the approximation in Eq.(37) 
and the mentioned ambiguity. After the integration over z in Eq.(34) we have neglected the iO+. 
Keeping this term the result reads: 

Tn^{2^^)gl^J'^'^^e-'^-'^^^-'^^ (53) 

With this term we can expand the matrix element in A and perform the integration over A ana- 
lytically. The result is: 

T ~ 2 g V 1 (//0; + A)|O,„.+,|/,0;)) 

~ ^Mme (^c - i0+)2"+l (p+)2n+2 ' ^^^^ 

where operators Oj with J = 2, 4, 6, • • • are the standard twist-2 operators in the light-cone gauge: 



= J = 2,4,6,--., (55) 

and the J-th moment oj of gluon distribution function is then given by: 

{h{v)\Oj\h{v))^{v^yaj. (56) 

In the following discussion we neglect the spin-flip terms, these terms can be discussed in the same 
way and they do not contribute at the leading order. 

In the limit s ^ 4m^ ^ t we can first neglect the t-dependence of matrix elements and only 
keep the dependence of p"*" and of A"*". For A+ the matrix elements can be expanded as: 

(Mp + A)|02n+2|MP)) _^ +vr-^^^6 -a +Vx^fo r57) 

) k=i k=i 

using these expansions we obtain 

S 



(58) 



If we assume that for n — > cxd the behavior of 62n+2,fc is similar as a2„+2 or h2n+2,k converges 
faster than a2n+2i then the dominant contribution for small x^. comes from the first sum, i.e., the 
dominant contribution comes from the terms with the moments aj of the gluon distribution, and 
other sums with 62n+2,fc are suppressed by positive powers of Xc- The first sum can be written 
as a integration from, which is just the approximation used in Eq.(37). Now we consider the 
approximation in Eq.(39). Using translational covariance the integral can be written: 

Tn^{2m)gl^J^ -e-2^-=^'^-°^^(/i(p + A)|G"'+^(0)G'^';(-An)|/i(p)). (59) 

Similarly, the integral can be written as a sum: 

2 s ^ 1 (/^(p + A)|g"'+^(^a+)-G%|/^(p)) 
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We divide the sum into a parts with even n and another part with odd n: 

2 s ^ 1 + A)|G"'+^(^^+)^^^G'"'+^|/l(p)) 

4me ^0 (2a;, - i0+)2'=+2 (p+)2fe+3 " ^ ^ 

For A"*" ^ the matrix elements behave hke: 

(/i(p + A)|G°>+^(za+)2^G°.+^|fe(p)) _ 

(^Tj^^T^ = -(2fc+l)— a2.+, + 0((— )) 

= {2k + l)xca2k+2 + 0{xl) (62) 

Prom these one would neglect the matrix elements with odd numbers of derivatives, i.e., one 
would neglect the summation in the second line of Eq.(61), because the matrix elements with odd 
numbers of derivatives go to zero with A"*" — > 0, or with — > by noting A+ = —2xcP'^. Keeping 
only the leading terms of the matrix elements with even numbers of derivatives and neglecting the 
matrix elements with odd numbers of derivatives, one obtains: 

s °° 1 

Writing the sum into a integration form as done for Eq.(58), one obtains the approximation in 
Eq.(39) and we will get the result for Tr which is related to gh{x) with x = 2x,. = 4m,f./s. But, the 
matrix elements with odd numbers of derivatives in the second line in Eq.(61) can not be neglected. 
Although they are proportional to — > 0, but the power of Xc in the denominator in the 
second sum is 2/c + 2 in comparison with the power 2k + 1 in the sum of the first line. Keeping the 
leading terms for matrix elements with even numbers and odd numbers of derivatives, one obtains 

2 s ^ 1 2 s ^ {2k + l) 

Clearly both sums give equally important contributions to Tji. Actually every term in the expan- 
sion in A+ for matrix elements with odd numbers of derivatives can lead to a equally important 
contribution to Tr. To sum the leading order contributions one can use 

{h{p + A)\G'^'+^{d+rG'^'+^\h{p)) = ± ^,(^'':fe), (^)'(MP + A)|G"-+^ {Id+y G'^'+^lhip)) 

(63) 

to re-arrange the sum. With the same assumption made after Eq.(58) one will get the same result 
as that obtained with the approximation in Eq.(37). 

In this section we neglect corrections suppressed by the inverse of s in the limit s — > oo and use 
two assumptions to derive our results in Eq.(51) and Eq.(52): One is used in Eq.(37) and specified 
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in detail between Eq.(58) and Eq.(59), another is to use the assumed asymptotic behavior to 
determine the real part of the S-matrix element. Without these two assumptions one can relate 
Ffi{z) in the limit to the skewed gluon distribution by neglecting corrections suppressed by the 
inverse of s. The relation reads: 



Fr{z) ^ -g'^- — XiX2/g(xi,X2), (64) 
where we used the definition of skewed parton distributions in 0]. There exists different definitions 



and the relation between them is well discussed in [|T4|. The variables Xi and X2 are related to z 
as: 

xi = -{l-z), X2 = -{l + z). 65 

s s 

With this we end up with the result for Tr: 



with 



dX2J —^XiX2fg{Xi,X2), (66) 

[Xi — Z0+)(X2 + ^0+) 



Xi-X2 = -. (67) 

s 



Without any knowledge about fg{xi,X2), the integral can not be performed and numerical pre- 
dictions can not be made. 



With our results we are in position to compare our result with others obtained in [|TT[]. The 
starting point there is to consider the process in Eq.(l), where the initial photon is virtual. The 
interaction between the initial hadron h and the cc pair is taken as that between h and a cc dipole 
with a small transverse-size. Using the perturbative result, in which only the exchange of two 
gluons are considered, one obtains at the leading order of Q~^: 

^l*-^o ~ ^1(1 + ^Ps) N.(x)] p {1 + OiQ-')} , with x = ^. (68) 

To derive the result one uses a light-cone wave function to describe the nonperturbative properties 
of J/ ip and J/ ip should be taken as massless to consistently perform the expansion in 

In the 

above equation the terms relevant for our comparison are given explicitly. As also pointed in the 
leading log approximations in ln(l/x) and in \n{Q'^ / Aq^j^) are required to identify the quantity 
X = Q'^/s as the variable of the gluon distribution gh{,x). The produced J/ifj is longitudinally 
polarized. It should be noted that by taking the cc pair as a small dipole it is equivalent to neglect 
higher orders in (5^^[|ll|]. To include the production of a transversally polarized J/^/', the above 
result is re-derived by keeping the mass Mjj^ and it is found that at leading order of the 
result is obtained by replacing with + Mj^^. The results for the production of J/ijj with 
a real photon are obtained simply by taking ^ in the above result, after the replacement 



of with + Mjf^i, |1TT|, |T^. Several modifications of the results are then introduced, and the 



relation between the forward S-matrix element and the gluon distribution becomes complicated. 
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We will compare the result in [Tl|, |I2| without these modifications. The result in |Tl|, |12| by 
taking ^ is: 

dcT 1 Atti^ 

-^\t-.o ~ (4;;^l(l + '^«) [^9hix)] |2, with X = (69) 

where the term with jSs is the real part of the amplitude and is given by: 

^Ps[x9,{x)]^^l^^. (70) 

If we take a in Eq.(42) as a small parameter, this part is the same as ours. The main difference 
between our result and that given above is that the S-matrix element is related to the gluon 
distribution gh{x) at different x. It is possible that the limit — >■ can not be taken because 
higher orders in are neglected in Eq.(68). 

Similar derivation of the results by using a nonrelativistic wave-function for J/ip is also done 
T5[| , by keeping the mass of J/ip at the leading order of and neglecting higher orders in 



m 



Q . Setting Q = one obtains: 



da 16na^rjJll^e^ ^ (71) 

dt 3 aemM°ji^ 

This result looks similar as ours, but the differential cross section is related to xgh{x) at a; = 2xc 
as the result in Eq.(69) and the real part of the amplitude is neglected. 



4. Comparison with experiment 

In this section we will compare our results with experiment performed at HERA. In the last 
two sections we have given our results, but not specified the energy scale /i, at which the nonper- 
turbative quantities like the gluon distribution are defined. To compare with experiment one must 
choose this scale. In general one may take the scale as a soft scale at order of Aqcd because of 
the emission of soft gluons. If one takes as a soft scale, then the scale should be the same in the 
production of T. However, there can be exchanges of hard gluons between quark lines, between 
gluon lines and between quark- and gluon lines. The effect can be studied with perturbative QCD 
and log terms like ln(yu/mc) appear, one may then identify the scale fi as rric for J/ip and as for 
T, respectively. We will compare our results with experiment with these two choices. It should 
be pointed out that significant corrections to our results for J/ip exist, as discussed below. Before 
these corrections are under control, a detailed comparison of theoretical results with experiment 
can not be made. In this section we only give numerical results based on our results at leading 
orders without any modifications, suggested by those corrections. 

To predict the total cross section we assume that 

^W^^«^0)e-l': (72, 
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where the slope parameter B is measured as i? ~ 4.5GeV^^||T6|. If the scale /i is taken as a soft 
scale, there is no detailed information available for the gluon distribution at the scale, one may 
know that xgh{x) is divergent like when x 0. Then the cross section behaves like: 

a{J/^) = A[^^j , l^o = lGeV, (73) 

where W^p = y/s and 5 = Aa. If we assume that the parameter B is the same for T, then we can 
predict the cross section of T by determining the parameters in Eq.(73) from experimental data 
for J/-?/^!!^. This can also be considered as we neglect the /^-dependence. We fit the published 
HERA data with Eq.(73) and obtain: 

5 = 0.83 ±0.13, A = 1.37±0.84nb, /d.o.f. = 0.31. (74) 

The fitting quality is quite good indicated by the small value of x^/d.o.f.. The fitting cure with 
experimental data is also shown in Fig. H. With the determined parameters we can predict: 

a(T) ^ 2.7 X 10-2 nb, for W^p = 120 GeV, 

(t(T) ^ 3.2 X 10^2 nb, for W^p = 143 GeV. (75) 

However these predicted values are too small in comparison with the central values of experimental 
results shown in Fig. ^. Different reasons are responsible for this discrepancy. Firstly we note 
that there are large errors in experimental data, even large errors in determining the energy W^p. 
Secondly, there are several uncertainties in our predictions, e.g., relativistic corrections arising 
in the expansion in Eq.(12), corrections from higher orders in the expansion in m~^, and also 
the leading order determination of the matrix element through the leptonic decay in Eq.(17), 
etc.. Among them the effects from relativistic corrections and from the determination of the 
matrix element can be most significant. Because the c- or c quark moves in the J/tp rest frame 
with the velocity Vc, which is not very small, indicated by ~ 0.3, the relativistic correction 
can be significant. In our approach this correction can be systematically added and its study 
is underway. In the determination of the matrix element we have used the leading order result 
for the leptonic decay in Eq.(17). It is well known that the corrections from higher orders in 
are large |p!7|, Also, the relativistic correction to the result is significant . It is clear that 
without these corrections an accurate prediction can not be made for the production of J/ip- For 
the production of T these corrections are expected to be small. 

If we take the scale fi as the heavy quark mass, we can predict cross sections without any 
input parameter, because the gluon distribution at the scale is determined. We use the recently 
determined sets of parton distributions GRV98||, CTEQ5|ll, MRST99[||] and MRST98LO||]. 



We do not use the newest MRSTOl distributions for J/^jj, because it has an unphysical zero at 
H = rric in the gluon distribution and for small x the distribution turns to be negative p3]. Instead 



we use an old set of parton distributions!^. The predictions with experimental data|jT6| are also 



shown in Fig. 2 with leading order gluon distribution, and in Fig.3 with the next-to-leading order 
gluon distribution. With these distributions we can determine the contribution of the real part 
of the S-matrix element to the cross section. The contribution is at 12%-, 20%- and 26% level 
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Figure 2: The cross section of elastic photoproduction of J/ip versus W^p, the photon-proton 



center-of-mass energy. The data points are pubhshed HERA resuhs[16|. The full solid line repre- 
sents a fit of the form a{J/ip) oc W^p. The theoretical predictions of present work using various 
parameterizations of the leading gluon density in the proton at ttIc scale are also shown. 



with the leading order gluon distribution of MRST98LO, CTEQ5L and GRV98LO, respectively, 
while with the next-to-leading order gluon distribution of MRST99, CTEQF3 and GRV98 it is 
at 6%-, 8%- and 12% level, respectively. From these figures we can see that all distributions 
give the results roughly with the same feature that the cross section increases with the energy, 
numerical results with the next-to-leading distribution of MRST99 are close to the experimental 
results for large s, while the prediction with other gluon distributions gives a rather large cross 
section. Although a rather good description of experimental data for J/ip can be given with the 
next-to-leading order gluon distribution of MRST99, one should keep in mind that our results can 
have significant corrections as discussed before and that it is not consistent to use next-to-leading 
order gluon distributions with our results at leading order of as- We note that the cross section 
does not depend on the renormalization scale /x. Since we work only at a^-the leading order, the /x- 
dependence can appear at the next-to-leading order of i.e., at af. To use next-to-leading order 
gluon distributions, the corrections from the next-to-leading order of as should be also included 
and the /^-dependence at is eliminated. These corrections come from exchanges of hard gluons 
between quarks and gluons in Fig.l. and they can be large because of that at /i = rric is not a 
very small parameter. Unfortunately, these corrections are unknown. It is clear that a conclusive 
comparison without those corrections, i.e., the relativistic correction and one-loop correction, can 
not be made with experiment. 

For T our predictions are obtained by taking /i = rub = 5GeV. The discussed corrections 
should be small for the case with T. We indeed find that our prediction with leading order 
gluon distributions is in agreement with experiment, although the experimental errors are large. 
Our numerical results are shown with experimental dataPS in Fig.4 with leading order gluon 
distributions, and in Fig. 5 with the next-to- leading order gluon distributions. The contribution 
of the real part of the amplitude to the cross section is at 36%-, 43%- and 47% level with the 
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Figure 3: The cross section of elastic photoproduction of J/ip versus W-yp, the photon-proton 



center-of-mass energy. The data points are pubhshed HERA resuhs[16|. The full solid line repre- 
sents a fit of the form ai^J/ip) oc W^p. The theoretical predictions of present work using various 
parameterizations of the next-to-leading gluon density in the proton at rric scale are also shown. 



leading order gluon distribution of MRST98LO, CTEQ5L and GRV98LO, respectively. With the 
next-to-leading order gluon distribution of MRSTOl, MRST99, CTEQF3 and GRV98 it is at 20%-, 
23%-, 29%- and 33% level, respectively. For T productions another possible correction besides 
those discussed before is that ai fi = rub the finite skewness, which we have neglected to relate 



the S-matrix element to the gluon distribution, can lead a significant effect. Studies in p6|, 27 



show that this is the case. Instead using gluon distribution one may need to use the skewed gluon 
distribution to make numerical predictions. However, this distribution is not well known and 
numerical predictions cannot be made. It deserves a further study of corrections to our results 
both for J/ip and for T. 



5. Summary 



In this work we have studied diffractive photoproduction of J/ip, where we have taken the 
gluons exchanged between the initial hadron and the cc pair as soft gluons. We then have used 
an expansion in to study the effect of the exchange of soft gluons. Our results have been 
derived with the assumption the exchange of two gluons in a special gauge. But, they can also 
be derived without the assumption in an arbitrary gauge. This can be done by formulating the 
problem as the splitting of the photon into a cc pair in a background field of gluons, which varies 
slowly in the space-time. Our results for the S-matrix element consist of a NRQCD matrix element 
representing the nonperturbative effect related to J/ and a matrix element of two gluonic field 
strength operators which are separated in the moving direction of J/ip in the space-time. The 
matrix element of the field strength operators characterize the nonperturbative effect related to the 
initial hadron. In the limit of s — oo the forward S-matrix element is related to the skewed gluon 
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Figure 4: The cross section of elastic photoproduction of T(IS') versus W^p, the photon-proton 
center-of-mass energy. The data points are pubhshed HERA resuhs[|^. The theoretical predic- 
tions of present work using various parameterizations of the gluon density in the proton at rrib 
scale are also shown. 



distribution. With a reasonable assumption the forward S-matrix element can also be related to 
the usual gluon distribution, this enables us to make numerical predictions and to compare with 
experiment. 

Since we take the exchanged gluons as soft gluons, it is not very clear how to identify the 
renormalization scale n to make numerical predictions. However, with the knowledge about the 
asymptotic behavior of the gluon distribution gh{x) with x — > our results show that the total cross 
section increases with increasing energy. We first take yU as a soft scale and fit experimental results 
for J/ip with our results. The experimental data show clearly that the cross-section increases with 
the energy as W^. With determined parameters we can predict the cross section for T. However, 
the predicted cross sections are too small. Possible reasons are large corrections to the leading 
order results for J/ip. We also take /i as a hard scale, i.e., as the mass of the heavy quark, with an 
existing gluon distribution we are able to give a rather good description of describe experimental 
data for J/ip, although large corrections are neglected. These corrections for T are expected to be 
smaller those for J/ip. We indeed find that our predictions with leading order gluon distributions 
for T production are in agreement with experiment, but there are only two data points with large 
errors. It deserves a further study of corrections to our results and of experiment. 

Our approach is distinctly different than previous approaches. In previous approaches one 
start to use perturbative QCD to study diffractive production with a virtual photon in the initial 
state, where the virtuality of the photon is large to ensure that perturbative QCD can be used 
for gluon exchange. Keeping leading order in and the finite mass of quarkonia one obtains 
results for the cross section. Then the cross section of diffractive production of a real photon is 
obtained by setting = 0. Clearly, this setting can not be done properly, because higher orders 
in are neglected. In our approach we take the exchanged gluons as soft gluons and use an 
expansion in the inverse of the heavy quark mass to handle the exchange of soft gluons. Although 
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Figure 5: The cross section of elastic photoproduction of T(IS') versus W^p, the photon-proton 
center-of-mass energy. The data points are pubhshed HERA resuhs[|^. The theoretical predic- 
tions of present work using various parameterizations of the gluon density in the proton at rrib 
scale are also shown. 



the forward S-matrix element can be related to the gluon distribution as that in the previous 
approaches, but the relation is different, as discussed in detail in Sec. 3. 

In our approach corrections can be systematically added. Among them relativistic correc- 
tion can be most important for production of J/ip, and a study of the relativistic correction is 
underway for getting accurate results of theory. With the accurate results of theory diffractive 
photoproduction may provide an interesting way to study the nonpertubative nature of nucleon 
and of nuclei. 
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Appendix 



In this appendix we study the integral which appears in Eq.(48). We define 

/(/?)=/ rfAA-^(cos(A) +isin(A)) = / d\\-^e'^, 
Jo Jo 

where (3 is reaL We make an exchange of variable A^~^ = z, then the integral becomes 

I roo ^ 

I{(3 = 1 — a) = — / (i2;exp(i2;<i ). 
a Jo 



(76) 



(77) 




Figure 6: The integration contour in complex 2;— plane. 

We extend the integrand as a complex function and consider the contour integral where the contour 
in the complex 2;-plan is specified as in Fig. ^ with 9 = ^an. Using the contour integral we have: 



11/"°° i 

I{P = l — a) = exp{i-a7i) / dRexp{—R'^) 

a 2 Jo 

for a > the integral is finite and it is real. With this result we obtain 

/■oo ]^ roo 

/ ciAA~^sin(A) = tan(-a7r) / ciAA~^cos(A) 

which is used in Eq.(48). 



(7J 



(79) 
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